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Abstract 

Within the Minimal Supersymmetric Standard Model (MSSM) we systematically investigate the 
bounds on the mass of the lightest neutralino. We allow for non-universal gaugino masses and 
thus even consider massless neutralinos, while assuming in general that i?-parity is conserved. Our 
main focus is on laboratory constraints. We consider collider data, precision observables, and also 
rare meson decays to very light neutralinos. We then discuss the astrophysical and cosmological 
implications. We find that a massless neutralino is allowed by all existing experimental data and 
astrophysical and cosmological observations. 
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1 Introduction 



1.1 Motivation 

The LHC is scheduled to start taking data in 2009, experimentahy opening the window to the TeV 
energy scale [1]. One of the main goals is to search for physics beyond the standard model (SM). 
A promising candidate is weak-scale supersymmetry (SUSY), which among other attractive features 
stabilizes the hierarchy between the weak scale and the Planck scale [2]. Supersymmetry predicts a 
large number of new particles, which if kinematically accessible, should be observable at the LHC [3-6]. 

The lightest supersymmetric particle, the LSP, plays a special role in the search for supersymmetry 
at colliders. If weak-scale supersymmetry is realised in nature, the production rate of the heavier 
squarks and gluinos will be dominant, due to their strong interactions. But these heavier particles 
rapidly cascade decay to the LSP in the detector, typically with no detached vertex. Since this occurs 
in nearly all supersymmetric events, the nature of the LSP and its behaviour are decisive for all 
supersymmetric signatures at the LHC. Several LSP candidates have been discussed in the literature, 
e.g. the lightest neutralino, the gluino [7,8], the lightest stau [9-11], and the lightest sneutrino [12,13]. 
In the case of conserved proton hexality, Pg [14, 15] (or conserved -R-parity^ [16]) the LSP (and the 
proton) is stable. A stable light gluino LSP has been excluded by the LEP data [17]. A stable (left- 
handed) sneutrino LSP is also experimentally excluded [18]. Furthermore cosmologically, a stable LSP 
must be electrically and colour neutral [19]. This leaves as the most widely studied candidate the 
lightest neutralino: Xi! which is a very promising dark matter candidate [19-21]. Further possibilities 
beyond the MSSM are the gravitino [22] or the axino [23], which we do not discuss here. We mention 
as an aside that in the case of ensuring the proton stability via baryon-triality [24, 25] the LSP decays, 
and in principle any supersymmetric particle can be the LSP. When embedding such models in minimal 
supergravity the stau, the right-handed smuon, the right-handed selectron, the sneutrino, the lightest 
sbottom or stop, as well as the right-handed sstrange or sdown have been found to be consistent with 
experiments [26, 27] . 

In this paper, we focus on a particular aspect of the neutralino LSP, Xi? namely 'How light can 
it be?'. Our main interest hereby are the laboratory constraints, i.e. 'What can one learn about the 
neutralino mass from collider or fixed-target experiments7\ Thus we initially put aside all aspects of 
the neutralino as a potential dark matter candidate. We shall discuss the cosmological implications of a 
very light neutralino towards the end of this paper. One should keep in mind that potential cosmological 
constraints or requirements can be avoided by adding a small (or large) amount of i?-parity violation. 
Potential other dark matter candidates in this case have been discussed in the literature, see for example 
Refs. [28-30]. Unless otherwise indicated, we assume in our analyses that i?-parity is conserved. 

The particle data group (PDG) cites as the laboratory bound on the lightest neutralino mass [31] 

mj^o>46GeV (1) 

at 95% C.L., which is based on the searches at LEP (the limit quoted above was obtained in the 
analysis of the DELPHI collaboration [32]). This has been obtained by taking into account the LEP 
searches for charginos, which yield lower limits on the SU(2) gaugino mass parameter, M2, and the 
higgsino mass parameter, /x. Furthermore, this bound assumes an underlying supersymmetric grand 
unified theory, based on a simple Lie group. In this case, the U(l)y gaugino mass parameter Mi is 
related to M2 by 

Ml = ^ tan2 0^M2 « ^Ma , (2) 

so that the experimental bound on M2 implies a lower bound on Mi. Thus the chargino searches yield 
lower limits on M2 and /i and indirectly on Mi , which together give rise to a lower bound on the mass 
of the lightest neutralino, see Eq. (1), via the neutralino mass matrix, see Eq. (8) below. 

^This is equivalent for the renormalisable superpotential. 
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It is the purpose of this paper to investigate the consequences of dropping the assumption Eq. (2), 
which invahdates the bound in Eq. (1). Such a scenario could occur, for instance, in the case of an 
underlying string unification with a semi-simple gauged Lie group [33]. In this more general scenario, 
Ml and M2 are both free parameters. As we shall see, this additional freedom allows the LEP bound, 
Eq. (1), to be avoided. This raises the question of the corresponding new lower mass bound on the neu- 
tralino. Such models have been considered in Refs. [34-41], however while still requiring the neutralino 
to provide the entire dark matter of the universe. As stated above, our main focus here are the labora- 
tory bounds. We systematically demonstrate that a massless neutralino is consistent with all present 
laboratory data. We shall consider bounds from direct collider searches, precision electroweak data, 
and rare meson decays. We then consider the astrophysical (supernova) and cosmological implications 
of a light neutralino, with a particular focus on a massless neutralino. 

1.2 Outline and Connection to Previous Work 

In early supersymmetric model building, the superpartner of the photon, the photino, 7, was often 
considered to be massless or very light (mass below 1 keV) [42,43]. The vanishing mass was guaranteed 
by a global ii-symmetry [43], under which the parameters of the supersymmetry transformation, the 
Grassman variables 9, 9, transform non-trivially. However, such an i?-symmetry has several problems. 
First, it prohibits all gaugino masses, including the gluino mass. As mentioned above, a light gluino 
has been excluded by LEP data [17]. Second, spontaneous breaking of the global i?-symmetry leads 
to a problematic axion. The axion and the light gluino can only be evaded by large explicit breaking 
terms which however render the original symmetry meaningless [44]. It is possible to gauge a U(l) R- 
symmetry in N=l local supersymmetry [45]. However in order to cancel the mixed triangle anomalies, 
we must extend the field content by non-singlet fields under the SM gauge symmetry, or consider 
a family dependent U(l)ij [44], see also Ref. [46]. In contrast, we shall here consider a very light 
or massless neutralino, where the small mass comes about through a modest amount of fine-tuning 
between the parameters Mi and M2 and consider no further symmetry. This is discussed in detail in 
Sect. 2. It is then a question of phenomenology, which we will discuss in detail in this paper, to see 
whether such a light neutralino is consistent with all laboratory, astrophysical and cosmological data. 
We would also like to point out that moels with a light neutralino have recently been obtained in the 
context of gauge mediated supersymmetry breaking, where heavy messenger masses are provided by 
the adjoint Higgs field of an underlying SU{5) grand unified theory [47]. For related work on light 
neutralinos in other extensions of the MSSM see, for example, Refs. [48,49]. 

Following up on the early supersymmetric model building, the phenomenological implications of 
a stable light photino were mainly discussed in the context of cosmology, namely as a dark matter 
candidate^ [19-21]. Such a photino constitutes hot dark matter. If this would provide the entire dark 
matter, it is inconsistent with structure formation in the early universe [53-57]. More recently it was 
found [34-38,40,41] that a light cold dark matter neutralino of mass 0{5 GeV) can be obtained within 
the MSSM without Eq. (2), while still providing dark matter to be consistent with the WMAP obser- 
vations [58,59]. As opposed to the earlier work on a light photino, we consider here a light neutralino; 
in fact we show that these are predominantly bino, given experimental constraints. Furthermore we 
consider significantly lighter, possibly massless, neutralinos. A stable and light neutralino constitutes 
hot dark matter. For m^o <C 25 eV it contributes only a small amount to the energy density of 
the universe and should thus be consistent with observations. We determine the precise number in 
Sect. 6.2.1. Another dark matter candidate is then however required [60]. A light photino or neutralino 
of mass 0(10 MeV) could also be produced in supernova explosions [61-65].'^ We reinvestigate these 
astrophysical and cosmological questions particularly for a massless neutralino in Sect. 6. 

As mentioned above, the main focus of this paper is to analyse in detail the issue of laboratory 
constraints on a light neutralino. We consider collider searches for light neutralinos, the contribution 

■^See also the later work [50-52], where photinos of mass 0(1 GeV) were considered together with a light gluino. 
''For a particular application of light but massive photinos see also Refs. [66]. 



of light neutralinos to precision electroweak observables via radiative corrections and neutralinos in 
rare meson decays. At colliders, the most promising processes for direct searches are 

e++e- ^ xlxh, (3) 
e+ + e- ^ e+e-xlxi- (4) 

These were investigated early for a photino [67,68]. Three of the present authors [69,70] have recently 
presented the first complete calculation of the process Eq. (3) for a general neutralino, focusing on 
heavy neutralinos, and analysed the resulting phenomenology at the ILC including beam polarisation. 
In Sect. 3.1, we reanalyse this for the case of a massless neutralino. We also consider the case of 
associated production of the lightest and next-to-lightest neutralino [71]. 

A second method to investigate a light neutralino in the laboratory is via its virtual corrections 
to precision observables. A first point that must be considered when discussing light neutralinos is 
whether the light mass is stable under radiative corrections. The neutralino also contributes to all 
electroweak observables via radiative corrections. For early work see for example Refs. [72-74], and 
Ref. [75] for a recent review. In most analyses the assumption Eq. (2) was typically made; the case 
of varying the parameter Mi independently was not systematically analysed. In general one would 
expect that a very light neutralino could give rise to significant supersymmetric corrections to the 
precision observables. We study this issue in Sect. 4, where we discuss in detail the total and invisible 
Z widths, Tz and Tinv, the W boson mass Mw, the effective leptonic weak mixing angle sin^^effi the 
anomalous magnetic moment of the muon {g — 2)^ as well as electric dipole moments in the case of 
complex parameters. 

Concerning rare meson decays to neutralinos, the first process studied in the literature was [76-78] 

K^nx'ixl (5) 

For a light neutralino the experimental signature is equivalent to the SM decay K — > ttvl' , which 
has since been observed experimentally [79-82]. The branching ratio still has a large experimental 
error but is so far consistent with the SM prediction [83-86]. This can thus be used to set bounds 
on the related supersymmetric parameters. However, typically the decay involves virtual squarks and 
if these are sufficiently heavy, no bound on the neutralino mass is obtained. For a massless or near 
massless neutralino, we can consider this instead as a lower bound on the relevant sfermion mass, as 
in the original study of the anomalous magnetic moment of the muon [87]. In Sect 5, we present a 
detailed overview of i^- meson, pion, S-meson, and quarkonium decays [3,77,88-93]. A more complete 
treatment of the bounds resulting from meson decays, including complete (higher order) calculations 
of the decay rates is deferred to a separate paper [94] . 



2 The Neutralino Framework 



2.1 Tree- level and higher-order corrections 

We shall work in the minimal supersymmetric SM (MSSM). The part of the Lagrangian which describes 
the neutralino mixing is given by [5] 
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and the fermionic fields are two-component Weyl spinors [95]. Here 

/Ml -Mzs^ cos P 

M2 Mzc^ cos (3 

-Mzs^ cos (3 Mzc^ cos j3 

\ Mzs^sinP -Mzc^smf3 -fj, 



(6) 
(7) 



Mzs^ sin/3\ 
-Mzc^sva.13 

y 



(8) 



and 



i^l = (B, W", h\, hi) . (9) 



Mz is the Z gauge-boson mass and = sin^^, c„ = cos 0w = Myy/Mz, where is the weak mixing 
angle, /i is the higgsino mass parameter and tan/3 = ^2/^1 the ratio of the two vacuum expectation 
values of the Higgs doublets. 

The chargino mixing is described by the following matrix [5]: 

^ = -{tl;-fXi)+, where i;+ = {W+,hlf, = {W- Jjf , (10) 

/ M2 V2Mwsmp\ 
X = [ . (11) 

\V2 Mw cos p 11 ) 

The mass matrix X is diagonalised by a biunitary transformation: 

dmg{m ±,m ±) = U*XV-\ (12) 

with U , V unitary 2x2 matrices, see for example [3]. The lower experimental bound on the lightest 
chargino mass is [31] 

m-±>94GeV. (13) 

Al 

Scanning over the parameter space, taking this bound into account, leads to a lower bound on both 
1^1 and M2 [48] 

M2 > lOOGeV. (14) 
If the GUT relation between Mi and M2, Eq. (2), is assumed, then this implies 

Ml > 50 GeV . (15) 

The usually quoted lower bounds on the lightest neutralino mass arise from imposing the experimental 
bound on the lightest chargino mass, Eq. (13), in conjunction with the assumption of the GUT relation 
between Mi and M2, Eq. (2). The value given in Eq. (1) furthermore takes into account results from 
other searches for supersymmetric particles and constraints from the Higgs sector, see Ref. [32] for 
details. 

If instead the theoretical assumption of Eq. (2) is dropped and Mi and M2 are treated as inde- 
pendent free parameters, there is an additional freedom in determining the lightest neutralino mass. 
In fact we can require a vanishing lightest neutralino mass at tree-level, by setting the determinant of 
the mass matrix, Eq. (8), to zero 

det(A^^o) = 0. (16) 

This is equivalent to [96, 97] 

[M2MI4 sin(2/3) + Ml ( - M2/i + M|c^ sin(2/3)) ] = . (17) 
The solution /i = is excluded by the above chargino bounds. Solving for Mi yields 

_ M2M|sin(2/3)4 
' /iM2 -M| sin(2/3)4- ^ 

Accordingly, for fixed values of /i, M2 and tan (3 one can always find a value of Mi such that the lightest 
neutralino becomes massless. Typically, the second term in the denominator is much smaller than the 
first, so that the resulting expression for Mi is approximately given by [making use of sin(2/3) = 
2 tan (3 /{I + tan^ (3) w 2/ tan /?, for tan /3 > 3] 

^ Ml.n.msi ^ / JO /ISOGeV-N 

fi \tan p J \ /X / 
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Figure 1: Bino admixture of Xi (l^ft plot) and masses of charginos and neutralinos (rigiit plot) for M2 = 
200 GeV, tan/3 = 10, and Mi as given in Eq. (18), such that m^o = [98]. Left of the vertical lines at 
fi « 120 GeV, the chargino mass is m-± < 94 GeV. In the right panel, the dotted line indicates the kinematic 
reach of LEP2 {^/s = 208 GeV) for e+e" — s- XiXi production {i = 2,3,4), and the dashed line indicates the 
mass of the Z boson, Mz « 91 GeV. Note that X2 is nearly mass degenerate with Xi for n > 120 GeV. 

This typically implies Mi <C M2, fJ.. In this parameter region the lightest neutralino Xi is predominantly 
bino, i.e. it couples to hypercharge. The bino admixture is typically above 90% in the parameter 
range where the chargino mass bound is satisfied. The masses of the other neutralinos and charginos 
are of the order of M2 and fi, see Fig. 1. 

The results given in Eqs. (18) and (19) have been obtained at the tree-level. Since the chargino/neu- 
tralino sector is described by the three independent parameters Mi, M2 and /i (in the case of real 
parameters, and we furthermore assume that tan/3 is determined via the Higgs sector of the MSSM), 
only three of the six chargino and neutralino masses are independent. Consequently, the other three 
masses are predicted. When including radiative corrections, one has to choose a certain renormalization 
scheme to define the physical meaning of the parameters. The three independent parameters of the 
chargino/neutralino mass matrices can be traded for three masses that are specified as input quantities. 
While in general loop corrections can give rise to a shift between the physical mass and the tree-level 
mass, the three masses chosen as independent input parameters do not receive higher-order corrections 
by construction. For the discussion of a very light neutralino it is thus convenient to choose m^o 
as one of the input parameters in order to avoid that higher-order corrections could drive it away 
from zero. Such a scheme where the two chargino masses and the lightest neutralino mass have been 
renormalized on-shell has been applied for the calculation of higher-order corrections in the MSSM 
chargino/neutralino sector [99,100]. Accordingly, once the mass of the lightest neutralino has been 
arranged to be small at tree-level, with an appropriate choice of renormalization scheme it remains 
small also if higher-order corrections are taken into account. The same also holds in the case of complex 
parameters which is discussed in the next subsection. 

2.2 Complex parameters 

The condition for a massless neutralino can also be obtained for a CP-violating neutralino sector. Then 
the parameters 

Ml = [Mile^"^! and = [/ile^"^" (20) 

are complex and have CP-violating phases (pi and <j)^. We choose the convention where M2 is real and 
positive, absorbing its possible phase by a redefinition of the gaugino fields. In the presence of complex 



phases, two equations have to be separately fulfihed to have a zero mass neutrahno, 



Jm{det{M^o)} = , and 9\e{det{M^o)} = . (21) 
These conditions are equivalent to 

M|4sin(2/3)sin</.i-|//|M2sin((^i + </.^) = 0, (22) 

M2M|4sin(2/3) + |Mi|[-M2^fcos(0i + ,/.^) + M|c^sin(2/?)cos0i] = 0, (23) 

respectively, and they can be solved for the absolute values of 

I I M|c^ sin(2/3) sin a# . 2 /> sin((/>i + (/>^) 

W = X. ■ — — and Ml = -M2tan — — . (24) 

M2Sin(0i+0^) smc/)^ 

The equations can also be converted to any other set of two neutralino parameters, for example 

^ M|cism(2/?)sin ^i M|4 sin(2;3) sm ^i _ 

sm(0i + 0^) |/i|sm0^ 

Note that in the presence of non-zero complex phases, the conditions for a massless neutralino cannot 
always be fulfilled. For example, it follows from Eq. (25) that the phases have to fulfill sin^i/ sin((/)i + 
0^) > and sin / sin 0^ < 0. In the CP-conserving limit with vanishing phases, we retrieve the 
condition Eq. (18) for M\ for a massless neutralino from Eq. (23), with cos((/)i + 0^) — > 1, cos(^i — > 1, 
and Eq. (22) is trivially fulfilled. 



3 Collider Bounds 

In this section, we consider the bounds on a light neutralino from collider searches. We focus our 
discussion in particular on the direct searches performed by the experimental collaborations at LEP. 
The bounds from LEP are in general more stringent than the ones from other lepton colliders with 
lower energies and/or luminosities. We briefly discuss limits from 6-factories in Sect. 3.2. Concerning 
bounds from hadron colliders, at the Tevatron the large QCD background limits the sensitivity in the 
search for direct production of a light or nearly massless neutralino [101]. At the LHC there could 
be better prospects for detecting effects of a very light neutralino in cascade decays of other (heavy) 
SUSY particles [102]. 



3.1 Neutralino production at LEP 

If we assume a massless neutralino Xi by choosing M\ as given in Eq. (18), the mass of the next- 
to-lightest neutralino X2 mainly determined by the values of [i and M2. In Fig. 2(a), we show 
contour lines of the second lightest neutralino mass.^ Qualitatively, the dependence is m-j^o ~ M2 
for M2 <C ^, and similarly m^o ~ [i for M2 ^ /i. This can be also observed in Fig. 1, where we 
show the dependence of the masses on 11. Thus for [i or M2 < 200 GeV, the associated production 
of neutralinos, e+e~ X1X2' ^1- ^^^S- 3' would be accessible at LEP up to the kinematical limit of 
^ = m^o = 208 GeV, if m^o = 0. 

In order to compare with the results of the LEP searches we make use of the model-independent up- 
per bounds on the topological neutralino production cross section obtained by the OPAL collaboration 
in the searches at LEP with = 208 GeV [103], 

a(e+e- ^ x\xl) x BR(x^ ^ Zx?) x BR{Z ^ qq). (26) 

''Note that here and for the following scenarios we choose an intermediate value of tan {3 = 10 and discuss fi > only, 
since the masses and cross sections change only slightly of the order of 10% if we take larger values of tan/3, and/or 
negative values of /i. 
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Figure 2: (a) Contour lines in the /i-M2 plane of the neutralino mass m^o. In the grey shaded area the 
chareino mass is to-± < 94 GeV. The dashed line indicates the kinematical limit m^o = -^/s = 208 GeV at LEP2. 

" Xl •'12 ^ 

Throughout we have chosen Afi such that to^o = 0. The lightest chargino is nearly mass degenerate with X2 ^'^^ 



M2 > 200 GeV and ^ > 125 GeV. (b) 95% confidence limit on the cross section a{e 
Zxi) at = 208 GeV (taken from Rcf. [103], Fig. 10). 




In this analysis a stable lightest neutralino is assumed. In Fig. 2(b) the observed limit at 95% confidence 
level on a{e^e~ — > X1X2) ^ BR(x2 ~^ ^Xi) is shown (taking into account that BR(Z — > qq) ~ 70%) in 
the m^o-m^o plane [103]. For m^o = (and m^o < 190 GeV) one can roughly read off the upper limit 

^(e+e- ^ xixl) X BR(x^ ^ Zx?) < 70 fb . (27) 

We analyze this bound assuming conservatively that BR(x2 ^Xi) = 1- general this branching 
ratio can be significantly smaller than 100% since other decay modes like X2 ~^ ^Xi ^.nd X2 ~^ e^e^ 
(see below) can be open. Imposing the bound a{e~^e~ — > X1X2) < 70fb significantly constrains the 
parameter space, since the typical neutralino production cross sections are of the order of 100 fb for light 
neutralino and selectron masses. In Fig. 4(a) we show contour lines of the cross section cr(e^e~ X1X2) 
in the //-M2 plane for tan/3 = 10 and degenerate selectron masses mg^ = rrigj^ = nie = 200 GeV. We 
observe that there is a large region in the 11-M2 plane where the cross section is higher than 70 fb. 
In order to fulfill the bound, the selectron masses have to be sufficiently heavy. It should be noted 
that the bino-like Xi couples preferably to the br, which is exchanged in the t and u channels, see 
the Feynman diagrams in Fig. 3. Thus, the bound on the neutralino production cross section can be 
translated into lower bounds on the selectron mass = mg^ = mg, for m^o = 0. In Fig. 4(b), we 
show contours of the selectron mass, such that the bound a{e^e~ — > X1X2) < 70 fb is fulfilled. For 
example, for a fixed selectron mass of rrie = 200 GeV, the area below the 200 GeV contour in Fig. 4(b) 
has a cross section in excess of 70 fb. 
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Fi gure 4: (a) Contour lines in the plane of the neutralino production cross section uie^e — > X1X2) 

with tan (3 = 10, and mg^ = mg^ = = 200 GeV, at ^/s = 208 GeV. At each point, Mi is chosen such that 
™x? ^ Contour lines in the /i-M2 plane of the lower bounds on the selectron mass to^^ = m^^ = mg, 

such that cr(e+e^ X1X2) = 70 fb for m-^o = with tan/3 = 10. In (a), (b), the dashed lines indicate the 
kinematical limit m-^o = ^/s = 208 GeV, in the grey shaded areas the chargino mass is 'm^± < 94 GeV. Along 
the dot-dashed contour in (b) the relation mg — m^o holds. 



It should be noted that above the dot-dashed contour in Fig. 4(b), the selectron is lighter than 
the neutralino X2- Thus, in that region the two-body decay into a selectron and electron, X2 ~^ e^e^, 
is open, leading to a reduction of BR(x2 ^Xi)- The decay of the second-lightest neutralino into 
a selectron and electron is followed by — > e^Xi- The reduction of BR(x2 — ^ ^Xi) implies a 
decrease of the lower bound on the selectron mass in this region. This effect is of minor relevance, 
however, since the bound on the selectron mass in this region is already close to the absolute lower 
experimental bound on the selectron mass, mg ~ 75 GeV [31]. It is clear from our analysis that for a 
sufficiently heavy selectron mass a massless neutralino could not have been discovered at LEP. Already 
for mg = 200 GeV, we have a significant range of allowed parameter space, as one can see in Fig. 4(b). 

Finally we want to note that the bound on the neutralino production cross section that we have 
used is indeed rather conservative. Since the decay X3 — > Zxi would also lead to the final state 
e+e~ — > qq + the associated production of the next heavier neutralino, e^e~ X1X3) would 
increase the signal rate if m^o < y/s. Since we have approximately rn^a /i, see e.g. Fig 1, the 
reaction e^e~ — > X1X3 would have been kinematically accessible for /i < 200 GeV at LEP2 with a 
center-of-mass energy of ^/s = 208 GeV. Including the additional cross section from e'^e" X1X3 
production in our analysis, we find that for ji < 150 GeV the selectrons must now be heavier than 
1 TeV. For > 200 GeV there is no additional bound. 

3.2 Radiative neutralino production 

An additional search channel at LEP for a light neutralino would be radiative neutralino production, 
e+e~ XiXil- However due to the large background from radiative neutrino production, e~^e~ — > 
i^uj, we find that the (theoretical) significance^ is at best S ^ 0.1 for C = 100 pb~^ and ^/s = 
208 GeV [69,70]. In addition, cuts on the photon energy or angle do not enhance the significance, due 
to similar kinematic distributions of signal and background. We find a similar situation at 6-factories, 
which are operating at the 6-meson resonances, ^/s ~ 10 GeV. For example, with C = 100 fb~^ and 
mg = 200 GeV, we find a significance of 5 < 0.1, and signal-to-background ratios not larger than 
1%. Further, an identification of the signal 'photon plus missing energy' is difficult due to the large 

^We define the significance S as tlie number of signal events over the square root of the number of background events. 



Q 



photonic background from the abundant hadronic processes at the 6-factories [104]. Note that at the 
ILC, however, radiative neutrahno production would be measurable, due to the option of polarised 
beams, which allows to reduce the background and enhance the signal at the same time [69-71, 105]. 



4 Precision observables 



Electroweak observables have in the last decades played a key role in constraining the free parameters 
of the SM and the MSSM (see e.g. Ref. [75] for a recent review in the context of the MSSM). In 
the following we study the impact of a light or massless neutralino on electroweak precision physics. 
Among the key observables in the electroweak sector are the mass of the W boson Mw, the effective 
leptonic weak mixing angle sin^ ^cflF (both discussed in Subsect. 4.2), the electric dipole moments of the 
electron, the neutron, and the mercury atom, and the anomalous magnetic moment of the muon (5— 2)^ 
(for the latter see Subsect. 4.3). The total Z boson decay width Tz and the invisible Z decay width 
Finv are potentially very sensitive to a massless neutralino, as in such scenarios the additional decay 
channel Z — > XiXi is kinematically allowed. The resulting tree-level constraints on a light neutralino 
were first investigated in Refs. [64,71]. In a first step, in Sect. 4.1 we reanalyse the impact of a massless 
neutralino on the total width and the invisible width of the Z boson, including full one-loop and leading 
higher-order corrections. 



4.1 Total Z Width Tz and Invisible Z Width r;„v 

A light neutralino with mass m^o < Mz / 2 can give contributions to the total and the invisible width 
of the Z boson, in addition to the decay channels into SM fermions 

Tz = r^^sM + r^? ' 

decay into SM fermions 

Tinv = Tinv^SM + T-j^O . (29) 

decay into neutrinos 

Potentially the additional contributions due to T^o can be large if the neutralino has a considerable 
non-bino like component, i.e. a sizable coupling to the Z boson. In Ref. [106], the processes Z — > XiXi 
and Z ^ ff have been calculated at 0{a) and supplemented with leading higher-order terms from the 
SM and the MSSM (see also the discussion below). The corresponding results, which are the currently 
most accurate MSSM predictions for these quantities, are used in the following to analyse the impact 
of a massless neutralino on the Z decay width. The experimental values for the total width and the 
invisible width of the Z boson are [31, 107] 

r^^'P = 2495.2 ± 2.3 MeV , (30) 
r.^P = 499.0 ± 1.5 MeV . (31) 

Below, we label the experimental errors of these two quantities as cjp^ and a'^^ , respectively. In our 
numerical analysis, we show the results for 

^Tinv ^ Tinv - r;=^P , (32) 
STz^Tz-rf'', (33) 

i.e. the difference of the MSSM prediction and the experimental result for the invisible Z width, (^Tinv, 
and the total Z width, 6Tz. 

In the following, we investigate both (JFinv and 6Tz in two representative SUSY parameter regions. 
As a first scenario, we choose fairly light scalar fermions and set the diagonal soft SUSY-breaking 
parameter MsusY to 250 GeV (in this section we choose this parameter to be equal for all sfermions). 
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In Fig. 5, we show 5Tinv in the upper and dVz in the lower plot, as a function of M2 and fi. Mi is fixed^ 
via Eq. (18). The remaining SUSY parameters are tan /3 = 10, Ar = At = Ai, = nig = Ma = 500 GeV. 
Here Af (/ = t,b,T) denotes the trilinear couplings of the Higgses to the sfermions, nig is the gluino 
mass, and M^ denotes the mass of the CP-odd Higgs boson. The deviations from the experimental 
central values as given in Eqs. (32), (33), are indicated as experimental nx a contours of the respective 
observable. In addition, the 95%C.L. exclusion bounds of 'm~± > 94 GeV [31] on the chargino mass 
from direct searches are marked by dashed white lines. The allowed regions are towards larger values 
of IM2I and To the left, the plots stop at around fi w —230 GeV, where the square of the lighter 
scalar top mass turns negative. Fig. 5 clearly displays that for both observables the MSSM prediction 
can deviate considerably from the experimental values. This is in particular the case for small and 
small I M2 1 . Nearly all of the parameter space ruled out at the 5a level for Finv and at the 3a level for 
Tz is, however, already excluded due to direct chargino searches. For the interpretation of these plots 
it is furthermore important to keep in mind that the results for Finv and Tz do not only depend on 
fi and M2, but on all the other SUSY parameters as well. This means in particular that an apparent 
la effect can easily be caused or canceled out by, for instance, a change induced by MgusY, the mass 
scale of the sfermions, which is known to have a strong impact on the decay into SM fermions (see also 
the discussion in Ref. [106]). Furthermore even in the SM, Fjnv is predicted to be slightly larger than 
the experimentally measured value, resulting in a ~ lo" deviation. 

As a second parameter range, we consider heavier sfermions: MsusY = 600 GeV. The remaining 
parameters are given by tan /3 = 10, A^- = At = Af, = nig = Ma = 600 GeV. As before, we show in 
Fig. 6 the results for the invisible Z width Finv (upper plot) and the total Z width Tz (lower plot), as 
a function of M2. Mi is fixed in the same way as before. Here the plot extends beyond fi = —600 GeV. 
The theory predictions are again confronted with their experimental values as described above. One 
finds somewhat better agreement of the Tz predictions with the experimental data. This is mainly 
due to the higher mass scale of the scalar fermions, which leads to the prediction of lower Tz values 
for a considerable part of the SUSY parameter space (see Ref. [106]). Despite the presence of F^o, the 
overall prediction for Tz even reaches values below the experimental la range, as indicated by the Oct 
and —la contours. Finv is in slightly better agreement with the data than before, but a la deviation 
can be observed for the entire allowed ;U-M2-plane. 

In summary, Tz and Fjnv cannot exclude a massless neutralino. The parts of the //-M2 planes that 
lead to a large deviation from the experimental values are mostly already excluded by direct chargino 
searches. 



4.2 W Boson Mass M\;y, Effective Leptonic Mixing Angle sin^ 0eff and Tz 

Next we analyse the impact of light neutralinos on Mw and sin'^ ^eflf (and extend the Tz analysis). 
The most up-to-date calculations for these quantities were presented in Refs. [108] and [106]. These 
predictions contain the full 0{a) MSSM result, as well as all known MSSM corrections beyond one- 
loop order, in particular universal contributions entering via Ap of 0{aas) [109] and 0{afi^) [110,111]. 
Concerning sin^^eff, in order to reproduce the best available SM results [112,113] in the decoupling 
limit, also the full electroweak 0{a'^) SM results [114-118], mixed electroweak and QCD SM corrections 
ofO{aas) [119] and 0{aal) [120], as well as leading beyond two-loop SM terms [121-123] are accounted 
for. We show these results as a function of Mi in an interval covering the solution of Eq. (18), or 
equivalently given by 

M2M|4sin2/? 
Af2/x-Af|4sin2/3 

with 5 Ml ranging roughly from —100 GeV to +100 GeV. Appropriate choices of M2 and fj, allow to 
analyse a Xi with a substantial bino, zino or higgsino component. 



j"^^"frr^:r:. +^Mi, (34) 



'A negligibly small offset value is added to Mi to acquire numerical stability, while scanning the /i-M2-plane. 



1 n 



-200 200 400 600 

^ [GeV] 




-200 200 400 600 

^ [GeV] 



Figure 5: The difference of the experimental value and the theory prediction for the invisible Z width, 
(fFinv, (upper plot) and the total Z width, 5Tz, (lower plot) in the /z-M2-plane, both including the 
process Z XiXi- Deviations of the theory predictions from the experimental data are indicated as 
,5r\nv = (Finv - rf^P) = (10, 5, 3, 2, 1) X (Tif^ (upper plot) and dVz = {Tz - T^'^) = (20, 10, 3, 2, 1, 0) x 
(lower plot) contours. The SUSY parameters were fixed as tan/3 = 10, MsusY = 250 GeV, A,- = 
At = Ab = m-g = Ma = 500 GeV. For Mi we use Eq. (18) (see text). 



-600 -400 -200 200 400 600 

^ [GeV] 




-600 -400 -200 200 400 600 

^ [GeV] 



Figure 6: Difference of experimental value and theory prediction for invisible Z width, (JFinv, (upper 
plot) and total Z width, 6Tz, (lower plot) in the /i-M2-plane, both including the process Z — > XiXi- 
Deviations of the theory predictions from the experimental data are indicated as JFinv = (Fjnv — F°^^) = 
(10,5,3,2,1.5,1) X o-r^P (upper plot) and 5Fi„v = (Fz - F^^'P) = (20,10,3,2,1,0,-1) x a^""^ (lower 
plot) contours. The SUSY parameters are chosen as tan/3 = 10, MsusY = A-r = At = Ah = nig = 
Ma = 600 GeV. Mi is defined via Eq. (18) (see text). 



In a first scenario a bino-like x? is obtained by setting tan/3 = 10, MgusY = 250 GeV, A^- = 
At = Ab = fj, = rUg = Ma = 500 GeV, M2 = 200 GeV and Mi according to Eq. (34), with the 
above mentioned range for 5Mi. This results in m^o values between and ~100 GeV. What is of 
interest here is not the absolute value of Mw and sin^ O^s (and Tz), which again strongly depends on 
the remaining SUSY parameters, but the change induced by varying Mi and thus also m^o. Clearly 
also these effects are somewhat dependent on e.g. the sfermion mass scale, but the effects due to the 
neutralino sector are the dominant ones. The results are displayed in Fig. 7, showing the dependence 
of Mw and sin^ ^efr on Mi and m^o. For complementarity with Sect. 4.1, we have included also the 
result for Tz- We show the two cases with T{Z XiXi) either included or not included into the total 
Z width. The left column shows the results as a function of Mi , while the right column as a function of 
the (corresponding) values of fn^o. The positive and negative values for Mi result each in one branch in 
the plots on the right. As it happens for the chosen set of parameters, Tz and M\y lie fully within the 
experimental one sigma range, while sin^ 6es is outside the la boundaries everywhere in the plot. As 
mentioned above, the relative agreement between the MSSM predictions and the experimental results 
would have been modified if we had chosen different MgusY, At^b, or M2 (or other SUSY parameters). 
As can be seen in Fig. 7, none of the observables has a strong dependence on Mi. The shift induced 
in Mw, for example, is ~ 0.5 MeV for the full Mi range. Even with an anticipated ILC precision 
of 6M^^^ = 7 MeV [124, 125] this is a marginal effect. The situation is similar for sin^ O^s. Here 
one can observe a ~ 0.4 x 10~^ variation with Mi, which is again marginal even compared with the 
ILC accuracy of 5 sin 

(9ILC = 1,3 X 10-5 [126]. Also the effects induced in Tz, even when including 
the Z-decay to neutralinos, are insignificant compared to the experimental errors, due to the small 
coupling of the bino-like Xi to the Z boson. 

As a second representative scenario, we study the case where tan (3 = 10, MgusY = A^ = At = 
Ab = Ma = m-g = 600 GeV, = 125 GeV, and M2 = 200 GeV. Mi is set according to Eq. (34), 
where 5Mi = —100 GeV to +100 GeV. For this choice of parameters the positive and negative values 
for Ml result each in one branch in the plots on the right, ranging up to m^o < 70 — 90 GeV. The 
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smaller value of ^ here allows for bigger non-bino like Xi-components. The numerical effects induced in 
the precision observables are therefore expected to be bigger. This expectation is confirmed by Fig. 8, 
where as before the left column shows the results as a function of Mi, while the right column as a 
function of the (corresponding) values of fn^o. The variation in M\y is ~ 5 MeV, so roughly a factor 
ten bigger than before, sin^ ^eff shows a variation of ~ 2 x 10-^. Still, even in comparison with the 
anticipated ILC precisions these are barely observable effects. The largest effects in Fig. 8 can be 
observed for the prediction of Tz- The impact of the Z-boson decay into a pair of light neutralinos 
on the total Z-boson width is seen to be quite substantial for this choice of SUSY parameters. With 
today's experimental accuracy, a 2a effect is observable in the presence of a massless neutralino. 

To conclude, even the most precise anticipated measurements of My/ and sin^ 9^^ are not able to 
constrain the mass of x?- The only possible constraints originate from the Z decay width, although 
even these are not very powerful for most of the parameter space, as they require a rather sizable 
coupling of Xi to the Z boson. Most of the parameter space which fulfills this requirement is already 
ruled out by direct chargino searches, as discussed in Sect. 4.1. The constraints are also easily satisfied 
by increasing the bino component of the light neutralino. 

4.3 Electric Dipole Moments and Anomalous Magnetic Moment of the Muon 

The effect of a small or vanishing mass of the lightest neutralino on electric dipole moments (EDMs) 
and on the anomalous magnetic moment of the muon, [g — 2)^, is shown in Fig. 9. For the EDMs 
the one- and two-loop formulas of Refs. [127-130] have been used, see Refs. [131,132] for reviews. 
The SUSY contributions to {g — 2)^ = 2a^ are based on the one- and two-loop formulas given in 
Refs. [133-136], see Ref. [137] for a recent review. These calculations have been performed using the 
computer code FeynHiggs [138-141]. 
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Figure 7: Total Z width Tz in- and excluding the process Z — > XiXii Mh/, and sin^^eff. The SUSY 
parameters are chosen as: tan (3 = 10, -/V/susY = 250 GeV, Ar = At = Ab = = nig = Ma = 500 GeV, 
M2 = 200 GeV. Ml is chosen according to Eq. (34), where 5Mi is varied from —100 to 100 GeV, and 
the two branches in the plots on the right result from negative and positive Mi. The green shading in 
the plots for Tz and Mw indicates that the whole area lies within the experimental la range. 
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Figure 8: Total Z width Tz in- and excluding the process Z — > XiXi: Mw, and sin^6'cfr. The SUSY 
parameters are chosen to be tan/3 = 10, -MgusY = A^- = At = Ab = Ma = rrig = 600 GeV, fi = 
125 GeV, M2 = 200 GeV. Mi is set according to Eq. (34), where 6M1 = -100 GeV to +100 GeV, 
and the two branches in the plots on the right result from negative and positive Mi. The green shading 
indicates the experimental la range of the three observables. 
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For the various EDM measurements upper limits exist as given in Tab. 1. For the anomalous 
magnetic moment of the muon a "discrepancy" of the experimental result from the SM prediction has 
been observed [142] 

flexp _ ^theo ^ (27.5 ± 8.4) X 10-^°, (35) 

equivalent to a 3.3 a effect.'^ While SUSY contributions could easily explain this "discrepancy", a 
massless neutralino could in principle lead to a too large contribution to a^. 
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Table 1: Present bounds for EDMs [146-148]. 
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Figure 9: The EDMs of the electron (upper left plot), the neutron (upper right) and mercury (lower 
left), as well as the shift in ((7 — 2)^ normalised to {g — 2)^{Mi = 500 GeV) (lower right plot), are shown 
as a function of Mi (see text). The other parameters are Mff± = MsusY = 500 GeV, Af = 500 GeV, 
M2 = // = 200 GeV. For the EDMs we have fixed tan/3 = 10 and set the phase of Mi to 4>i = 
7r/50, 7r/20, vr/lO. For the shift in {g — 2)^ we fixed the phases to be zero, but varied tan /3 = 5, 10, 40. 
m^o = is reached around Mi ~ 0. 



^Tliree other recent evaluations yield slightly different numbers [143-145], but similar discrepancies with the SM 
prediction. 
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The SUSY contribution to the EDMs and to {g — 2)^, shown in Fig. 9, have been evaluated as a 
function of Mi with the other parameters set to MsusY = — ^00 GeV, = 500 GeV, M2 = /i = 

200 GeV. For the EDMs we have fixed tan/? = 10 and set the phase of Mi to (pi = 7r/50, 7r/20, vr/lO. 
For {g — 2)^, we fixed the phases to be zero, but varied tan/3 = 5,10,40. According to Eqs. (18), 
(19), (24), (25) a massless neutrahno is reached around Mi ~ 0. We show the contribution relative to 
{g — 2)^(Mi = 500 GeV), which is 0(10"^) for this set of parameters. 

The SUSY contributions to the EDMS involving the lightest neutralino go to zero for vanishing 
m^o, as can be seen in the upper left, upper right and lower left plot in Fig. 9, where we show the EDM 
of the electron, the neutron and mercury, respectively. Consequently no lower bound on m^o can be 
set. For large values of cpi, on the other hand, an upper limit on Mi and thus on m^o can be derived 
as can be seen in the upper right and lower left plot of Fig. 9. The variation of (g — 2)^, shown in the 
lower right plot of Fig. 9, stays below ~ 0.5 x 10^^'^ and is thus well below the current uncertainty of 
8.4 X 10^^'^ [142]. Therefore, while (for this exemplary set of parameters) the SUSY contributions is 
of the same size as the deviation between the SM and the experimental result, no experimental limit 
on m^o can be set. 
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5 Rare Meson Decays 

In the SM, '^5i mesons (vector mesons) can decay into a neutrino pair, e.g. 

J/ii{lS),T{lS) ^ ViVi, i = e,//,r, (36) 

where J/ip and T denote cc and hh (ground) states, respectively. Below we also investigate (p = ss 
and the light mesons p and a;, which are superpositions of uu and dd states. An example Feynman 
graph for a meson decay into a neutrino pair at the parton level is shown in Fig 10(a). If we allow for 
neutrino masses, which enable a chirality flip, also ^5*0 mesons (pseudoscalar mesons) can decay into a 
neutrino pair, e.g. [149, 150] 

■K^ ^ UiVi i = e,ii,T , (37) 

where again Fig. 10(a) represents the tree-level graph. In both cases, the neutrino pair remains 
unobservable. It is possible to look for the invisible decay of mesons, Eqs. (36), (37), for example via 
the decay chain 

i){2S) JUj + 'K-n; i = e,p,T, (38) 

by tagging on the invariant mass of the di-pion system. To our knowledge, this idea was first proposed 
in Ref. [151], as a test of the SM neutral current, and has since been widely employed [152-155]. In 
Ref. [156] it was proposed as a method to look for physics beyond the SM, then a light gravitino. It 
has since been used to look for many different aspects of beyond the SM physics [89-93, 157-159]. 
Mesons can also decay to a lighter meson and a neutrino pair, e.g. [80,82] 

-fC"*" — > 7r"'"fjZ/i, i = e,ii,T. (39) 

A parton level Feynman graph is shown in Fig. 10(b). The rare decay, Eq. (39), has a similar signature 
to the cascade, Eq. (38), and has in fact been observed [80,82]. The partial decay width is consistent 
with the SM [160,161], albeit with a still large experimental error. It has also been used to investigate 
beyond the SM physics [76-78,162-165]. 

In the Pg conserving MSSM with a stable, light neutralino, we can replace the neutrino pair in the 
decays Eqs. (36), (37) and (39) by a pair of light neutralinos. The latter would also remain unobservable 
and we would obtain similar signatures, possibly slightly modified by the light neutralino mass. We thus 
have potential new decay modes which can be used to search for, or constrain, the light neutralino 
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(a) Meson decay into a neutrino pair. (b) decay into a neutrino pair and a pion 
Figure 10: Example diagrams for meson decays into a neutrino pair 

scenario [3,76,77,88-92,157,163]. The mass range of the decaying mesons {e.g. ^ 130 MeV, 
niK ~ 500 MeV, mx ~ lOGeV) enables a test of various light neutralino masses [76]. 

Motivated by our discussion in Sects. 2 and 4, we assume throughout this section that the neutralino 
is pure bino, denoted as B, and very light, i.e. that the neutralino decay mode is kinematically 
accessible to the meson. We recall that a pure bino has no interaction with the Z boson and its 
couplings to right-handed sfermions are stronger than to left-handed sfermions. The corresponding 
supersymmetric Feynman diagrams for the meson decays can be quite different than those for the SM 
decays. We list a few examples in Fig. 11. 




(a) pion decay (b) Bs decay (c) J/tp decay 

Figure 11: Example diagrams for the decay of a vr, Bs, and J/'i/'-meson into a bino pair. 

In the following, we give an overview over the present bounds existing in the literature, modifying 
them to binos, where necessary. The analyses carried out so far involve significant approximations. An 
up-to-date, detailed analysis of meson decays to light neutralinos, in particular including the necessary 
higher-order calculations, is beyond the scope of this paper and will be presented in Ref. [94]. However, 
given the current status of the data, cf. Table 2, we expect the present conclusions to be robust. 

5.1 Pseudoscalar Mesons 

Pseudoscalar mesons (psm), e.g. tt^, K, r], rj' , and Bg can in principle decay into a light neutralino pair. 
The related photino decay vr'^ — > 77 was first computed in Refs. [3,77]. Later, in Ref. [88], the decays 
7]-0 — > 77 and TT^ — > 777 were investigated. If we rescale the result in Ref. [88] from a squark mass 
of rriq = 70GeV to a squark mass of rriq = 300 GeV [166], we obtain an upper bound on the photino 
branching ratios, which is well below the current experimental limits for the neutrino decays [31] 

BR(7r° ^ z^P) < 2.7 • lO"'^ , (40) 
BR(7r° ^ Z/P7) < 6 • 10"^ . (41) 

This conclusion also holds for the corresponding bino decays. For their theoretical estimates the 
authors in Ref. [88] assumed a large left-right mixing in the squark sector, not taking into account 
flavour changing neutral current constraints. Due to the very small supersymmetric branching ratio, 
no bound is obtained on the bino mass. 
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In Ref. [3, 77] no left -right mixing was assumed. We reestimate their branching ratio for a pure 
bino and for heavier squarks and assume that the relevant squarks are degenerate. We also generalize 
to an arbitrary psm. The corresponding Feynman graph for a pion is shown in Fig. 11(a). The decay 
width to a bino pair is then given by 



r(psm ^ BB) = -^C(psm)^J^ Jl - ^ . (42) 
4cos^f7 " 



mj- V m? 



"q V ""psm 

The constant C depends on the considered meson and is given by 

-0\ _ 2 J^2 



C{7^n = 9/.0 , (43) 



C7(r/) = - (sin - ^ cos 9sh) , (44) 



C(r?') = - (sin + ^ cos Osfs) , (45) 
= ^16^ ' 

where f^o, fi, fs, and are the decay constant of the vr'^, r/i, 773, and Bg, respectively; 6*1^8 are the mixing 
angles between the rji and r^g states. Vcb is the charm-bottom CKM matrix element. These are all 
given in Refs. [31,167]. c parameterises the flavour mixing in the squark sector; representative values 
for c and /b^ are taken from Ref. [89]. 

At the parton level, the decay Bg — > BB proceeds via b — > sBB, cf. Fig. 11(b). The flavour 
changing neutral current (FCNC) is possible because the left-handed squark mass matrices and the 
quark mass matrices are not (necessarily) simultaneously diagonal [89] . This mismatch is parameterised 
by the constant c. The authors of Ref. [89] obtain for a 80 GeV squark mass a branching ratio of about 
O(10~^ — 10~^), depending on the amount of flavour violation. In Table 2, we have rescaled their 
result and present the branching ratio for a squark mass of 300 GeV. 

In Table 2, we list in the upper part the numerical values for the branching ratios for the various 
meson decays to binos. Here we have assumed rUg = 300 GeV. Note that the decay width, Eq. (42), 
peaks at a value = nipsm/V^- In order to obtain conservative bounds in Table 2, we assume in 
turn for each meson this value for the bino mass. It can be seen that the branching ratios for a meson 
decay into a light bino pair are several orders of magnitude lower than the experimental bounds on the 
respective invisible widths. This is mainly due to the small B mass in the numerator (because of the 
helicity flip required for a decay of a pseudoscalar into a pair of binos) and the large squark masses in 
the denominator of Eq. (42), suppressing this weak decay mode. Consequently no bounds on a light 
B can be set from these decay modes. 

5.2 Vector Mesons and Quarkonium Decay 

The decays of vector mesons to photinos were first considered in Refs. [3, 157, 169]. The decay of the 
various T(nS') (excited) states into a neutralino pair and its discovery possibility at B factories is 
discussed in detail in Ref. [92]. 

In Ref. [169] , the decay of quarkonium into gluinos and photinos has been calculated. Theses decays 
violate parity. The coupling strength of photinos and gluinos to left- and right-handed squarks are 
equal. Consequently, the masses of left-handed and right-handed squarks must be different to allow for 
these decays. However, the first argument does not hold for the decay into binos: The coupling strength 
to right-handed particles is larger than to left-handed particles, i.e. parity is explicitly violated. Hence, 
this decay is also possible when left- and right-handed squark masses are equal to each other. 

The quarkonium decay width can be approximated by, cf. Ref. [172], 
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Table 2: Comparison of maximum values of various branching ratios of meson decays to light binos. 
We show the theoretical maximum value and the current experimental bounds on the BR. For the 
theoretical computation we choose the bino mass to maximize the BR. The upper (middle) part shows 
the results for the pseudoscalar (vector) mesons, the lower part shows a rare decay. For the pseudoscalar 
mesons, the squark masses are set to 300 GeV. For the vector mesons, see Eqs. (48)~(52). All data 
are taken from Refs. [31,167] and from the specified literature. " — " indicates that no experimental 
bound has been derived yet. We have included also for comparison a related experimental result on 
the invisible decay of the T [92] . 



where my is the mass of the vector meson V = '^Si{qq), and rrig is the mass of the exchanged squark. 
Cv is a constant depending on the sum of the squared quark-squark-Bino-coupling strengths. For the 
different mesons, Cy is given by 
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Cp = j^Qggg , rriu^ = m^^ oo, rriu^ = mj^ = ruq = 300 GeV (48) 
25 

^ 72 ' "^"^ ^ ~^ "^"^ ^ = "i<? = 300 GeV (49) 
= rus^ oo, mg^ = rUg = 300 GeV (50) 

1 IT 

Cj/^ = — , ^ oo, mc-^ = mg = 300 GeV, = - (51) 

1 vr 
Ct = — , oo, = mg = 300 GeV, ^1=2' 

where we indicate the approximations we have made for the squark masses. In Eqs. (48)-(52), we 
have also listed the input parameters (squark masses and mixing angles) employed to derive the upper 
bounds on the branching ratio. They are chosen in such a way that the branching ratio is maximal for 
a massless neutralino. For the n, d, and s squark, we assume that left- and right-handed squarks do 
not mix and that all squark masses are degenerate. For the c and b squark, we include squark mixing. 
6c and 0^ denote the mixing angles in the scalar charm and bottom sector, respectively. 

As in the pseudoscalar meson case, the decays are strongly suppressed by the large squark masses. 
The results of the theory evaluation are compared to the experimental data in the middle part of 
Table 2. It can be seen that no bounds on the mass of a light bino can be deduced. The discrepancy 
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Figure 12: Examples for loop induced supersymmetric Kaon decays 

between the experimental bound and the theoretical estimate is quite large. We thus expect these 
results to be quite robust. More detailed formulae can be found in Ref. [94]. 

5.3 Loop-induced Meson Decays to Binos 

In Refs. [76-78,163], the related photino decay — > 7r"^77 was analysed. It can proceed at tree-level 
via the cascade decay 

E:+ ^ 7r+7r° (53) 
^ 77 

It can also proceed directly at the loop level, for which example Feynman graphs are shown in 
Figs. 12(a)-12(c). It was found that for large parts of the MSSM parameter space this decay is 
suppressed relative to the SM decay — > tt'^uu. Since the SM event rate is barely observable with 
present experiments, no bound on the photino mass or the relevant sfermion masses in the propagators 
is obtained. However, these results were derived when the correct mass of the top quark was unknown. 
Furthermore, only an incomplete set of the relevant one-loop Feynman graphs was evaluated, and the 
neutralino was restricted to be a photino. Nonetheless, the estimate obtained using the published 
results is well below the experimental bound, and we do not expect any bound from present data for 
a bino. However, in Ref. [94] this issue will be investigated based on a complete analysis. 

A similar analysis can be carried out for the decay B~^ — > ttBB. The corresponding SUSY box 
diagrams involve bosons, squarks, charginos, and Higgs bosons. As evaluated in Ref. [76,77], the 
suppression by the scalar fermions in the loops (which are assumed to exceed the experimental bounds 
of ~ 300 GeV) leads to a contribution that is too small to set any limits on a light bino. 

Another rare decay in which a light bino could in principle play a role is 6 ^ 57. However, at the 
one-loop level neutralinos do not contribute. The effect of a light B at the two-loop level would be 
well within the current theoretical and experimental uncertainties [173]. 



6 Astrophysical and Cosmological Bounds on the Neutralino Mass 

In this section we briefly consider the implications of a very light neutralino for astrophysics and 
cosmology. We focus on supernova cooling and the dark matter of the universe. In Ref. [174] the 
implications for a moderately light neutralino with m^o > 5 GeV for Big Bang nucleosynthesis are 

A 1 

discussed, which we shall not further consider here. 
6.1 Supernova Cooling 

During a supernova explosion neutrinos are abundantly produced in the dense core. They diffuse out 
with a time scale of 0(10 sec) [175-177]. The neutrinos with this time structure have indeed been 
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observed after Supernova 1987a [178, 179]. If light neutralinos exist with mass less or of order the 
supernova core temperature, Tc = 0(30 MeV), they can also be produced abundantly during core 
collapse. Depending on their interactions, these neutralinos escape freely from the supernova, rapidly 
cooling the core [180]. As the temperature drops, the neutrino scattering cross section drops with the 
square of the temperature, leading eventually to free-streaming neutrinos. Thus rapid cooling of the 
supernova with a time scale well below 10 sec is excluded by the neutrino observation from Supernova 
1987a. This can be used to set restrictions on the light neutralino mass, as well as its interactions. 
This was originally addressed for photinos in Refs. [61-63]. The authors of Ref. [64] derived important 
lower bounds on a light bino, which we briefly summarize here; see also Ref. [65]. Our focus here is on 
a massless neutralino. 

The two main neutralino production mechanisms in a supernova are electron-positron annihilation 
and nucleon-nucleon (NN) "neutralino-strahlung" : 



Once produced, the neutralinos have a mean-free-path, A^o, in the supernova core which is determined 
via the cross sections for the processes 



as well as the electron and nucleon densities. If A-j^o is of order of the core size, Rc = 0(10 km), or larger, 
the neutralinos escape freely and thus cool the supernova rapidly. However, if the neutralinos have 
masses m^o much greater than the supernova core temperature Tc, then their production is Boltzmann- 
factor suppressed and they affect the cooling negligibly, independent of A-r^o. Demanding that m^o be 
large enough that neutralino-cooling does not markedly alter the neutrino signal — particularly its time- 
structure — allows one to set a lower limit on the neutralino mass. Note that this limit depends strongly 
on the squark and selectron masses, which enter to the fourth power in Eqs. (54), (55) and Eqs. (56), 
(57) through the relevant propagators. 

A proper treatment of this problem would be to expand the existing supernova code(s) to include 
the production and the scattering of neutralinos. Thus the neutralinos would be involved in the com- 
plete time evolution of the supernova, which could affect the particle densities within the supernova 
and the supernova temperature as a function of time. Such a treatment is beyond the scope of this 
paper. A good estimate of the effect of the neutralinos on the supernova evolution can be obtained, 
if we use the existing codes [181-184] and treat this non-supersymmetric supernova as a fixed back- 
ground, i.e. we assume the neutralino effect on the evolution to be small. Using the resulting electron 
and nucleon densities, we can compute the production and scattering of the light neutralinos. We then 
employ the Raffelt criterion [185] , requiring that the maximal emitted energy from the supernova via 
neutralino radiation is < 10^^ erg. In Ref. [64], it was then found that for selectron masses in the range 
300 GeV < rUe ^ 900 GeV neutralino masses below 100 MeV are excluded. As the selectron mass is 
increased from 900 GeV the lower bound an the neutralino mass gradually decreases. For selectron 
masses above 1.2 TeV there is no lower bound on the lightest neutralino mass. 

Similar, however much less restrictive arguments also hold for the squark mass dependence. For a 
massless neutralino, squark masses between 300 GeV and 360 GeV are excluded. 

For selectron and squark masses below 300 GeV the mean-free-path of the neutralino is smaller 
than the supernova core size: A^o < Rc, i.e. the neutralinos are trapped, and diffuse out, just like the 
neutrinos. In this case the above approximate procedure is no longer valid and the neutralinos must be 
included in the numerical supernova simulation. This has to-date not been performed. Thus at present, 
massless neutralinos are not excluded by the Supernova 1987a observations for mg, ruq < 300 GeV or 
both me > 1200 GeV and rriq > 360 GeV. 
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It should be pointed out that to-date it is not yet possible to successfully simulate a full supernova 
explosion, in particular the outgoing shock wave still stalls. Thus some ingredient is still missing. 
More recently the simulations are being extended to three dimensions with the inclusion of turbulent 
effects [186,187]. An eventual full solution could in principle lead to a modification of the above results. 



6.2 Dark Matter 

In the MSSM with conserved proton hexality or conserved i?-parity, the lightest neutralino is stable 
and will contribute to the dark matter in the universe [19, 188-191]. If it is very light, i.e. relativistic 
at freeze-out, it will contribute hot dark matter. If the dark matter candidate is non-relativistic at 
freeze-out, it contributes to the cold dark matter of the universe. The structure formation in the early 
universe is best described by cold dark matter alone [192]. Thus the contribution of hot dark matter 
to the energy density of the universe is severely restricted. We discuss the resulting bound on the light 
neutralino mass, the Cowsik-McClelland bound, in Sect. 6.2.1. 

For cold dark matter, it is well known that for smaller interaction cross sections the resulting 
relic density is larger. Furthermore the reaction rate also decreases with decreasing mass of the dark 
matter candidate. Therefore, if we assumes that the candidate particle provides the required cold dark 
matter of the universe, we obtain a lower bound on the particle mass, the Lee-Weinberg bound. We 
shall discuss the resulting bound for a neutralino in Sect. 6.2.2. This bound assumes the standard 
big-bang cosmology with a radiation dominated universe prior to nucleosynthesis. If one drops this 
assumption, it has recently been shown that the lower bound is substantially weakened [41]. If for 
example one generalizes the MSSM to the next-to- minimal supersymmetric Standard Model (NMSSM), 
which contains an extra singlet chiral superfield, the bounds can also be weakened [48,49,92]. We 
restrict ourselves to the MSSM and the standard big-bang cosmology. 



6.2.1 The Cowsik McClelland Bound 

Here, we consider the case of a (nearly) massless neutralino, m^o < 0(1 eV). As for meson decays 
above, we restrict ourselves to a pure bino neutralino [71]. Since the very light bino contributes to the 
hot dark matter of the universe, we assume here implicitly that the cold dark matter originates from 
another source, see e.g. Ref. [60] for a review of alternate candidates. For a sufficiently light bino the 
contribution of the hot dark matter to the energy density of the universe is expected to be consistent 
with present observations. We wish here to quantify this statement and thus determine an upper mass 
bound on a very light stable bino. 

Note that if we were to assume that proton hexality or i?-parity is violated, then for such a light 
bino [m^o < 0(1 eV)], kinematically the only open decay mode \s B ^ (provided a lighter neutrino 
exists), which proceeds via a one- loop diagram. Using the computation in Ref. [193], we can estimate 
the lifetime in this case as 

,-,._Li„-»,ecf^;;Lyfi^y. (58) 
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where A is the relevant i?-parity violating coupling [26] , and m j is the mass of the sfermion entering 
the loop. For A = 0.01 (a typical upper bound), mj = 100 GeV and a bino mass = 1 eV, we obtain 
a lifetime of about 10^^ sec, well above the age of the universe. 

For a stable bino, the bino relic energy density, pj^, divided by the critical energy density of the 
universe, pc, is given by [194] 



Gn and Hq denote Newton's gravitation constant and the present value of the Hubble constant, 
respectively. Q{2>) is the Riemann zeta function evaluated at 3. is the photon temperature. Recall 



the effective internal degrees of freedom of a particle are given in terms of the internal degrees of 
freedom, g, by 




5efr = <! 3 ^ ^ . . (60) 
jg, for a lermion 

For the bino, which is a Majorana fermion, we have: gef[{B) = 2 x (3/4) = 1.5. 5*5 is given by 

9*s = 9i- +\ 5* • (j) ■ (61) 

i=bosons j=fcrmions 

Here, gi is the number of internal degrees of freedom and Tj the temperature of the particle species i, 
respectively. The sum runs over the index i for all species in thermal equilibrium at temperature T. 

In order for the bino hot dark matter not to disturb the structure formation, we conservatively 
assume its contribution to the total energy density of the universe to be less than the upper bound on 
the energy density of the neutrinos, as determined by the WMAP data [58, 59] 

Vt^h^ < [17^/i^]max = 0.0076 . (62) 

Light binos decouple at T ~ 0(1 — 10 MeV). This temperature is somewhat higher than the 
temperature where the neutrinos decouple, because the selectron mass is larger than M^y. The higher 
selectron mass leads to a smaller bino scattering cross section, i.e. earlier freeze-out. However, the 
temperature is well below the muon mass, so it is not necessary to know the exact value. Nevertheless 
at this temperature, we have 2 bosonic and 12 fermionic relativistic degrees of freedom (one photon, 
one Dirac electron, three left handed neutrino species, and one light Majorana neutralino) leading to 
g-tS = 12.5. From Eqs. (59) and (62), we find the conservative upper bound 

< 0.7 eV. (63) 

Thus a very light bino with mass below about 1 eV is consistent with structure formation. This line 
of argument was originally used by Gershtein and Zel'dovich [195] and Cowsik and McClelland [196] 
to derive a neutrino upper mass bound, by requiring Q.^ < \. We have here obtained an upper mass 
bound for a hot dark matter bino. 



6.2.2 The Lee Weinberg bound 

In this subsection, we determine the lower mass bound on a light neutralino from the requirement that 
it alone provides the required cold dark matter in the universe. This is based on the original work for 
massive neutrinos [197-200] and the resulting bound is referred to as the Lee-Weinberg bound in the 
literature. As a light neutralino, we consider here a neutralino with a mass below the LEP bound, 
Eq.(l), but which is nevertheless non-relativistic in the early universe, i.e. neutralinos with m^o > 
0(5 GeV). For simplicity, we shall again restrict ourselves to the case of a pure bino, as we are mainly 
interested in qualitative statements. This neutralino mass range has been investigated numerically in 
Refs. [34-40] in some detail. We present here a semi-analytical treatment of the Boltzmann equation 
following Ref. [194]. This gives some insight into the dependence of the neutralino relic density on the 
main parameters in our case: the bino mass ?n^, and the slepton mass m^. At the end of this section, 
we compare our results to Refs. [34-38,40]. 

For simplicity, we shall consider only the neutralino annihilation into leptons 

BB ^ il, I = e,iJL,T,i>e,Vn,VT ■ (64) 

The r-lepton is considered as massless, all sleptons are assumed to have common mass m^. The 
relevant annihilation cross sections are then related by 

a{BB £-1'^) = 16 a{BB r'tj^) = 16 (t{BB u^Vi) . (65) 
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Figure 13: Constant contours of the relic density of a bino type LSP as a function of the LSP mass, 
m^, and the slepton mass, m^. In the red (dark) shaded area the rehc density is in the ahowed 
range 0.091 < Q^Mh'^ < 0.129. The grey (hght) shaded areas are excluded due to LEP searches, 
< 80 GeV and the theoretical requirement: < m^. 




The cross section averaged over the thermal distributions of the incoming particles is given by [194,201] 
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Here T is the temperature of the universe. Implementing our specific cross section into the procedure 
outlined in Chapter 5.2 of Ref. [194], we obtain for the bino freeze-out temperature, Tf, 
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Here the thermally averaged cross section is evaluated at x = 1, i.e. T = m^. The effective number of 
relativistic degrees of freedom g^ is given by 
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(68) 



H{m^) is the Hubble parameter for T = m^. The bino contribution to the normalized energy density 
of the universe is then given by 
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(69) 



Here m-p\ denotes the Planck mass. We have plotted in Fig. 13 contours of constant relic density in 
the mg~m^-plane. The lower right-hand triangle of the figure is excluded since here the sleptons are 



lighter than the neutrahno, contrary to our assumption. The bottom horizontal band, < 80 GeV, 
is excluded due to LEP searches [31]. The red (dark shaded) band denotes the preferred relic density 
region from the 5-year WMAP data [59] 



n^h^ = 0.1099 ± 3an = 0.1099 ± 0.0186 , (70) 

where we have implemented the 3 sigma bounds, and where the 1 sigma error is an = 0.0062. 

As expected, we see in Eq. (69) that the resulting relic density is inversely proportional to the bino 
annihilation cross section. Furthermore 
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B 

Thus for fixed VLj^h? the slepton mass oc ^Jfn^. This can be seen in Fig. 13. Since observationally 
there exists an upper bound Q,j^h? < [rijj/i^]max = 0.129, this translates into a lower bound on the bino 
mass m^. We obtain the most conservative bound for the smallest allowed slepton mass = 80 GeV. 
Inserting these numbers we obtain 

m^>13GeV. (72) 
From Fig. 13 it should also be clear that one obtains an upper bound on for = and 

m^<419GeV. (73) 

This is then also the upper bound on the slepton mass in this scenario. 

The neutralino relic density has been considered extensively in the literature, see for example 
Refs. [19-21, 188-191]. More recently the question has been raised precisely how light the lightest 
neutralino can be, while still providing the entire (cold) dark matter in the universe as required by 
the WMAP data [34-39]. The authors allowed for non-universal gaugino masses, cf. Eq. (2), as 
well as a general neutralino admixture. They took into account all kinematically allowed annihilation 
products, as well as co-annihilation and resonant annihilation. In Ref. [34] a lower bound of about 
18 GeV was found for a non-relativistic neutralino. The lowest lower bound of about 6 GeV was 
found in Ref. [37,38] in models with a light pseudoscalar Higgs A with mass Ma < 200 GeV. In 
Refs. [36, 38] a semi-analytical approximation was also performed. They focused on the s-channel 
Higgs exchange specifically in the low region. It is encouraging how well these numbers agree with 
our approximate results, which only considered the neutralino annihilation to leptons. This is mainly 
because the resonant and co-annihilation effects we have ignored are not relevant in the low-mass 
region. 

In a further study, it was found that when allowing for explicit CP violation in the Higgs, and thus 
the neutralino sector the lower bound is reduced to 3 GeV [40]. 



7 Summary and Conclusion 

In this paper we have studied mass bounds on light neutralinos from collider physics, precision observ- 
ables, meson decays, and astrophysics and cosmology. We have focused on the question how light the 
Xi can be, in particular whether an essentially massless Xi is allowed by the experimental and obser- 
vational data. Assuming a semi-simple gauged Lie group the parameters Mi and M2 are independent. 
A very light neutralino is possible by mildly tuning Mi, resulting in a substantial bino component of 
the X?. 

Specifically, we have analysed the processes e"'"e~ — > X1X2 e+e" XiXi7 for a massless 
neutralino. The LEP searches have determined an upper bound on the cross section of the first process. 
This reaction depends on the selectron masses, due to the t-channel selectron exchange contribution. 
Consequently, we can translate the upper bound on the X1X2 production cross section into a lower mass 
bound on the selectron masses, depending on the choice of fi and M2. Assuming degenerate selectron 



masses and making the conservative assumption that the branching ratio of X2 i^^to x? a-nd a Z boson 
is 100%, for II < 150 GeV we find mg > 1 TeV. For fj, > 200 GeV no bounds on mg are obtained. As 
a consequence, this search channel yields no bounds on the mass of the lightest neutralino. Radiative 
neutralino production, e~^e~ — > XiXiT? which has been searched for at LEP2 and 6-factories, also does 
not result in any bound on m^o. 

Electroweak precision observables such as the W boson mass, M^y, the effective leptonic weak 
mixing angle, sin^ O^fi, or the invisible Z width, T-mv are potentially sensitive to a very light neutralino. 
For M]y and sin^^eff; where the neutralino contribution is a loop effect, we find that even the most 
precise anticipated measurements will not be able to constrain the mass of x?- The contribution to the 
invisible Z width, arising from the decay of the Z boson into a pair of very light neutralinos, can give 
rise to a measurable effect if the bino component of the Z boson is not too high. Also this constraint, 
however, is not very powerful for most of the parameter space, and the largest effects occur in parameter 
regions that are already excluded by the LEP searches for charginos. We have furthermore analyzed 
the effect of a light neutralino on electric dipole moments and the anomalous magnetic moment of the 
muon. Both do not show any relevant sensitivity to m^o. 

Next we investigated the bounds from rare meson decays. Pseudo-scalar mesons and vector mesons 
can decay into a pair of massless neutralinos. This decay, however, is suppressed by (large) squark 
masses and, in the case of pseudo-scalar mesons, a small neutralino mass. This decay has the same 
signature as the SM decay into a neutrino pair. Using the corresponding experimental upper bounds no 
constraints on m^o could be established. We also found no bound on m^o from the rare, loop-induced 

cascade decay — > vr'^XiXi- 

We next considered the astrophysical bounds on the neutralino mass from neutralino production in 
a supernova explosion. We found that a massless neutralino is consistent with the neutrino observations 
of supernova 1987a, provided the selectron mass (which enters in the production and scattering cross 
sections) is larger than 1.2 TeV and the squark mass is larger than 360 GeV. However, for squark and 
selectron masses below 300 GeV the neutralinos are trapped in the supernova, similar to the neutrinos. 
This scenario has not yet been sufficiently analyzed, so that also in this case a massless neutralino 
cannot be ruled out at present. 

Finally, we considered cosmological bounds on the neutralino mass. If the neutralino constitutes 
hot dark matter it is in agreement with the WMAP observations and structure formation for a mass 
below 1 eV. For the case that the lightest neutralino constitutes the entire cold dark matter of the 
universe, we performed a semi-analytical analysis and found a lower bound of m^o > 13 GeV. This 
is in remarkably good agreement with the more complete numerical computations which find a lower 
bound of about 6 GeV. 

Overall, we have found that a massless neutralino is consistent with all present laboratory exper- 
iments and astrophysical and cosmological observations. At the upcoming experiments at the LHC 
very light neutralinos could be produced in cascade decays of heavier supersymmetric particles. The 
anticipated precision in the determination of the differences between squared masses could allow to 
further test scenarios with very light neutralinos. Further insights can be expected from high-precision 
measurements at a future e'^e~ Linear Collider. 
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